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The 0(4) linear sigma model of the chiral transition in QCD is similar to models of the superfluid 
transition in “^He and ®He. Observations of vortex formation in superfluid helium have recently 
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held belief that this density is very small in comparison with the pion density. 
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I. INTRODUCTION 

Although QCD is believed to be the underlying theory 
of the quarks and gluons which constitute normal bary- 
onic matter, very little is currently known for certain 
about the nature of its chiral phase transition. A purely 
gluonic SU{3) theory is believed to have a first order 
transition Q|. In a world with just two massless quarks, 
however, or in other words an infinitely massive strange 
quark, the chiral transition is plausibly second order [Q. 
If the strange quark were massless the transition would 
again be first order Q . It is even conceivable that there 
is no real phase transition at all and that the change of 
state occurs by a smooth cross-over |^. For the physical 
values of the quark masses the nature of the transition is 
not certain although it is frequently assumed to be sec¬ 
ond order. It is also often assumed that the chiral and 
deconfining transitions are one and the same and are re¬ 
ferred to as simply the QCD transition. Here, we will 
adopt both these conventional assumptions and explore 
the consequences of the resulting similarity between the 
QCD transition and the superfluid phase transition. 

Collisions of highly relativistic nuclei offer the possibil¬ 
ity of producing quasi-macroscopic regions of dense nu¬ 
cleonic matter at a sufficiently high temperature that it 
might be possible to observe the QCD transition exper¬ 
imentally. In such ’heavy-ion collisions’, one typically 
collides Pb, Au, S or O at energies between 3 and 200 
GeV per nucleon. The high energies of the incident nuclei 
mean that, in the centre of mass frame, they are highly 
Lorentz contracted and resemble two pancakes approach¬ 
ing each other along the beam pipe at almost the speed 
of light. In fact, the energies of the incident nuclei are 
sufficiently high that the whole system is approximately 
Lorentz boost invariant. The practical result of this is 
that all physical quantities tend to be functions of proper 
time and hence observables are rapidity independent. 

Immediately after the collision, the two nuclei recede 
in opposite directions down the beam-pipe, leaving a re¬ 
gion of hot quark-gluon plasma between them. This then 
expands and cools through the QCD phase transition. 
Eventually the energy density becomes sufficiently low 


that the plasma hadronises to produce the pions, nucle¬ 
ons and kaons observed. Unsurprisingly by far the ma¬ 
jority of the products are pions since they are so much 
lighter than anything else. 

If one is to use such heavy-ion collisions to probe the 
nature of the QCD phase transition, it would clearly help 
to have some indication of what the experimental conse¬ 
quences of QCD ought to be. Due to the difficulty of 
applying conventional perturbative techniques to QCD 
or simulating the transition numerically, this is an in¬ 
credibly hard problem. Nonetheless there has been much 
work on possible experimental signatures, involving ob¬ 
servables such as the photon and dilepton fluxes and the 
K/tt ratio. There has also been interest in the effect 
of hydrodynamic instabilities during the cooling of the 
plasma. Possibly the most clear signature so far consid¬ 
ered, however, is the deficit of neutral pions which would 
arise from large regions of misaligned QCD vacuum act¬ 
ing as pion lasers, more usually known as dis-oriented 
chiral condensates. 

Another, very natural possibility, would be to look at 
the baryon density. Although the baryon number of the 
incident nuclei would probably almost all be contained 
in the receeding pancakes which constitute the remnants 
of the original particles, it is conceivable that there could 
be a significant baryon number density in the central re¬ 
gion immediately after the transition. Since the speed of 
the plasma is proportional to the distance from the col¬ 
lision point, this central region is also known as the cen¬ 
tral rapidity region. Any baryons in it would therefore 
have characteristically low longitudinal velocities com¬ 
pared with the receeding nuclei. 

Nonetheless, the central rapidity region is usually as¬ 
sumed to be baryon free, partly on the basis of string 
models, although there has until now been no work to 
predict the proton and neutron distributions directly [Q . 
It turns out, however, that in the context of the lin¬ 
ear sigma model the evolution of the baryon density has 
many similarities with the production of topological de¬ 
fects dl. Recent theoretical progress in this area [^,0, 
supported by experiments in superfluids [8[-p^, means 
that it is now possible to address the question of the 
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baryon density immediately after a heavy ion collision 
more directly. The resulting estimate of the initial baryon 
density is consistent with the conventional belief that the 
central rapidity region is almost baryon free. 


II. A MODEL FOR THE QCD TRANSITION 

In order to calculate the baryon density in the central 
rapidity region immediately after the transition, we must 
first choose a tractable model for its dynamics. We have 
already assumed the chiral and deconfining transitions to 
be one and the same and second order. Let us further 
assume that it is then reasonable to use two flavour QCD 
instead of the full theory. If we use superscripts L and R 
to distinguish the left and right handed sectors of QCD, 
then the breaking of the full chiral symmetry group:- 

U^{Nf)xU^{Nf) = Uv{l)xUA{l)xSU^{Nf)xSU^{Nf), 

where the vector C/v(l) and axial C/a( 1) sub-groups cor¬ 
respond to multiplying the left and right quark spinors 
by equal and opposite phases, to the residual symmetry:- 

SU{Nf)L+R 

would be described by non-vanishing expectation values 
for operators of the form:- 

M] = itLQRj)- 

Here q^i and qnj are left and right handed quark spinors. 

The question of what model to use to describe the 
dynamics of such an order parameter is far from clear 
cut. The full QCD lagrangian might be correct in princi¬ 
ple but makes calculation too hard since the interacting 
quanta are strongly interacting. Since below the critical 
temperature we have a good idea of what to expect phe¬ 
nomenologically, one often tries to deduce the more phys¬ 
ically relevant weakly-interacting degrees of freedom and 
use these to construct a more tractable model. A famil¬ 
iar example occurs in condensed matter systems where 
one starts with strongly interacting atoms in a crystal 
lattice and transforms this into a description in terms of 
weakly interacting phonons. For QCD, the weakly inter¬ 
acting perturbative degrees of freedom are usually taken 
to be mesonic, or in other words the pions. Exact trans¬ 
formations between the original quark-gluon degrees of 
freedom and the effective mesonic degrees of freedom are 
not known, but there are two phenomenological models 
commonly used, the Skyrme model 0 and the 0(4) lin¬ 
ear sigma model ||^. 

Both of these theories roughly reproduce multi-pion 
scattering amplitudes to order p^, or in other words at 
tree level. This is equivalent to treating the models as 
classical lagrangians. If one were to treat these new phe¬ 
nomenological theories classically, one would expect the 
baryons to be solitons. The justification for a classical 
treatment has been much discussed in the literature [Q. 


In fact, the classical approximation can be very similar 
to a fully quantum mechanical, although still approxi¬ 
mate treatment Physically, this is because the long 
wavelength modes of the scalar fields grow exponentially 
and their correlation length becomes larger than all other 
length-scales in the problem, including the inverse pion- 
masses and 1/Tc. In other words they become more 
classical as the transition progresses. Pictorially, one can 
think of classical, long wavelength ocean swells coming 
to dominate over the short wavelength chop as the tran¬ 
sition progresses. 

Let us first consider the Skyrme modeh- 

where the vector (p = {a, tt) is an 0(4) multiplet of real 
scalar fields, the vector tt representing the three pions 
and the a a sigma particle too massive to be observed at 
low energies. 

Due to the extra scale coming from the four deriva¬ 
tive term, this model has stable texture-like topological 
defects, usually called Skyrmions. It is possible to show 
that these textures have baryon number one and spin 
one half and they are therefore identified with the pro¬ 
tons and neutrons. Unfortunately, since exact isospin 
invariance was assumed there is no distinction between 
the two. Since experimentally it would be far easier to 
detect protons than neutrons this is potentially a seri¬ 
ous problem. The minimum energy solution for such a 
Skyrmion is of the form:- 

r 

TT = -/(r) sin0(r) 

a = f{r) cos0{r) (1) 

where / is some function of r which has to be calculated 
numerically by minimising the energy to give the result 
shown in figure ^ 

The conserved topological current associated with 
these textures is given by:- 

^ 2^4 4>a d^pb dxpc Oppd, 

where is the totally antisymmetric symbol in 

Minkowski space and Cabcd is the equivalent in 0(4) field 
space. The zeroth component of this gives the topolog¬ 
ical charge density, or in other words what turns out to 
be the baryon number density:- 



Since protons and neutrons are indistinguishable in this 
model this corresponds to the density of nucleons minus 
the density of anti-nucleons. 


-Cskyrme — 4’a+ ^^2^2 
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FIG. 1. This configuration of the 0(4) scalar field is stable 
for topological reasons. Such configurations are generically 
called topological defects and this particular example is called 
a Skyrmion. It corresponds to either a proton or neutron. 

From the point of view of predicting the baryon den¬ 
sity immediately after the QCD transition, however, this 
model has one fatal flaw, namely the Skyrme term itself. 
This term is clearly not conformally invariant since it is 
designed specifically to provide a scale for the protons 
and neutrons. Hence, a model incorporating such a term 
can not describe a renormalisation group fixed point such 
as a phase transition and in particular couldn’t describe 
the chiral transition. Another way of seeing the same 
thing is to regard the Skyrme term as a lagrange mul¬ 
tiplier which fixes the vacuum expectation value of the 
field. If this is constrained to be finite, clearly one can’t 
describe the symmetric phase in which, by definition, the 
vacuum expectation value should be zero. 

The only other alternative, without going to higher or¬ 
ders in some form of derivative expansion and ending up 
with a model which is totally impractical for calculation, 
is the 0(4) linear sigma model:- 

^ ^ - v‘^f + Ha. 


Again (f> = {a, tt) is an 0(4) multiplet of real scalar fields, 
TT representing the pions and the a a sigma particle too 
massive to be observed at low energies. A priori, this 
is a quantum field theory and the zero temperature val¬ 
ues of the parameters A, H and should be chosen to 
give reasonable agreement with experiment at low en¬ 
ergies. Given the other approximations associated with 
our choice of model, the exact experimental data used 
to choose the values of these parameters is not critical. 
Here the following values are chosen :- 

V = 87.4MeV 
H = (119MeV)3 
A = 20 

which are consistent with = 140MeV, rria = 600MeV 
and the pion decay constant /tt = 92.5MeV. There are 
other equally valid possibilities, however . With these 
parameters, the sigma model gives a reasonable descrip¬ 
tion of the phenomenology at less than IGeV. The phase 
transition takes place at roughly Tq ~ /^r and low en¬ 
ergy TT — TT scattering amplitudes come out about right. 
Probably the largest criticism of the model is that the 
a, which is far more massive than the quasi-goldstonian 
pions, has never been seen. 

One advantage of the linear sigma model is that all 
the critical exponents of the theory can be calculated, in 
both the static and dynamic renormalisation group, in 
the limit of small H. In fact, for the three dimensional 
theory, the linear sigma model is plausibly in the same 
universality class as the 0(4) Heisenberg ferromagnet ||] 
whose indices have previously been calculated to seven 
loops @]- Given the other approximations made here, 
this would seem to be plenty. In the conventional nota¬ 
tion, the critical indices are as follows:- 

a= 2-diy = -0.19 ±0.06 
P= ^{d-2 + r]) = 0.38 ±0.01 
7 = ( 2 - 77 )^ = 1.44 ±0.04 

HIZJI = 4.82 ± 0.05 

d — 2 ± 77 

z/= 0.73 ±0.02 (2) 

where d is the spatial dimension which we take to be 
three, notwithstanding the apparent flatness induced by 
the approximate Lorentz boost invariance, since the spa¬ 
tial structure of the field will later turn out to be impor¬ 
tant. If we define the relative temperature e = 1 — T/Tq 
in the conventional way, then the correlation length of 
the scalar field will be:- 



where ^0 ~ 0.7fm. 
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Similarly, the relaxation rate of the pion field may 
be obtained from dynamical renormalisation group 
arguments:- 

T = ToC 

where 


which unsurprisingly is of the order of the light crossing 
time of a pion. In three dimensions, the critical index 
z = d/2 = 3/2. 

The biggest draw-back with this model as far as we are 
concerned comes from Derrick’s theorem, which tells us 
that any renormalisable field theory involving only scalar 
fields can’t support stable solitons. In other words, the 
fact that there is no fourth order derivative term in the 
sigma model to provide a length scale for the baryons 
means that the configurations of the field which would 
be defects and represent the baryons are not topologi¬ 
cally stable. Whereas the Skyrme model had stable ob¬ 
jects corresponding to protons and neutrons but could 
not be a good description of the phase transition since 
it broke conformal invariance, the linear sigma model is 
conformally invariant at the transition but does not have 
stable protons and neutrons. 

The most obvious solution to this dilemma would be 
not to require the model to be a renormalisable field the¬ 
ory. Since we are looking for a phenomenological model 
rather than a fundamental field theory there is some jus¬ 
tification for this. In this case, however, it would not 
be clear exactly what was meant by conformal invari¬ 
ance and it would not be easy to calculate with whatever 
terms were necessary to ensure stable solitons. 

Another possibility would be to exploit the fact that 
we are only really interested in a finite sized volume of 
quark-gluon plasma and hence have boundary conditions 
so that Derrick’s theorem doesn’t necessarily apply. Cer¬ 
tainly, if we were to use the sigma model on a 2-sphere, for 
example by exploiting the axial symmetry to reduce the 
problem to just the radial and beam-pipe co-ordinates 
and then imposing the condition that the field should be 
zero at the edge of the plasma, there would be stable 
solitons, even though these would lack a scale and the 
nucleon size would not be fixed. 

Neither of these solutions is particularly appealing or 
clearly better than the other, however. Assuming that 
there really is a transition and that it is second order, 
then ideally one would like a renormalisable theory which 
is conformally invariant at the critical point and con¬ 
serves baryon number written in terms of the pion fields. 
It is not clear, however, that this is possible since the 
bosonic theory does not include either glueballs or all 
the quark flavours. It could also be the case that in real¬ 
ity there is a smooth cross-over rather than a real phase 
transition, in which case the theory wouldn’t have to be 
conformally invariant at the critical point. 


Here, however, we have assumed a second order tran¬ 
sition in two flavour QCD. In order to describe the tran¬ 
sition we therefore have to use the sigma model. In both 
the Skyrme and sigma model, however, the 0(4) scalar 
field represents the same physical degrees of freedom 
We know from the Skyrme model what configuration of 
the pion and sigma fields corresponds to a baryon and 
hence the same configuration ought also to correspond 
to a baryon in the sigma model. Indeed, if we con¬ 
sider the sigma model in thermal field theory, at low 
temperatures, the field would effectively be confined to 
the vacuum manifold, the texture configurations would 
be effectively stable and the model would be equivalent 
to the Skyrme model. We will therefore assume that 
the QCD phase transition can be described by the lin¬ 
ear sigma model and that the protons and neutrons are 
represented by Skyrmion-like configurations of the field 
notwithstanding the fact that they are not topologically 
stable. 

In fact, Skyrmion-like configurations of the scalar field 
will tend to collapse at the speed of light. If we were to 
take this seriously this would imply violation of baryon 
number. This is clearly a flaw in the sigma model and 
presumably arises since this model doesn’t contain all 
of the relevant physics. Is it, however, serious for the 
dynamics of the phase transition? Certainly, the pres¬ 
ence of topological defects can produce non-perturbative 
effects which can, for example, change the critical tem¬ 
perature 0- The time-scales for this process, however, 
are such that the rate of texture decay, or in other words 
baryon violation due to the inadequacy of the model, is 
always slower than the time-scale for breaking the sym¬ 
metry. The time-scale for the symmetry to be broken 
is of order 2Tl/Xv‘^c^ ~ Co/c. This should be compared 
with the minimum texture collapse time of ^/c. Since a 
texture corresponding to a proton or a neutron will al¬ 
ways be larger than the cold coherence length, it is safe 
to assume that textures will take longer to decay than 
the time available during the course of the transition. 
In fact, detailed studies of texture dynamics suggest that 
texture unwindings may be quite rare and even less likely 
than a simple time-scale argument suggests ||^. It will 
therefore be assumed that the of the 0(4) field is unlikely 
to be much influenced by the texture unwinding events 
which would describe proton / neutron decay. 

In conclusion, in order to calculate the baryon density 
immediately after a heavy ion collision, we will use the 
non-linear sigma model with critical exponents calculated 
in the H = Q limit to model the dynamics and assume 
that baryons correspond to configurations of the 0(4) 
scalar field which look like Skyrmions. 


‘Although in the case of the Skyrme model there is some 
ambiguity in the choice of which of the components is the 
massive sigma. 
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III. INITIAL CONDITIONS AND 
HYDRODYNAMICS 

Let us now consider the initial conditions. If the energy 
released in a heavy-ion collision is roughly equivalent to 
that produced in nucleon-nucleon collisions then one ex¬ 
pects the energy density immediately after the collision to 
be of order 3 GeV / fm^. Unlike many cosmological phase 
transitions, however, QCD is a strongly coupled system 
and, in the sigma model, the dimensionless coupling con¬ 
stant A must be of order 20 in order to approximate low 
energy pion cross-sections. At these sort of energy den¬ 
sities, this means that the interaction time-scale is likely 
to be far shorter than the cooling rate and the deconfined 
quark-gluon plasma is likely to come rapidly into thermal 
equilibrium at a temperature of a few hundred MeV. 

For example, if following Bjorken one assumes the 
initial energy density to be somewhere between 1 and 10 
GeV / fm^ and distributes this energy at about 400 MeV 
per quantum, then the mean free path works out to be 
roughly:- 


Am/p ~ f -] X (0.05 — 0.5)fm 

\ ^int / 


with a corresponding thermal equilibration time of the 
order of 1 fm / c. Thus, while the usual assumption of 
an initial thermal state is likely to be wrong in the com¬ 
paratively weakly coupled cosmological models, it may 
be a reasonable approximation for a heavy-ion collision. 

Typically people assume you reach 200 and 300 MeV at 
between 1 and 4 fm / c |15|. Here, as an initial condition, 
we will assume that the quark-gluon plasma is in thermal 
equilibrium at a temperature of 200-300 MeV at a proper 
time T/ « 1 fm / c after the collision. 

There is in fact another reason not to try to use 
the sigma model at times less than 1 fm / c. If we 
were to treat the sigma model as a phenomenological 
quantum theory rather than as a classical field theory, 
some momentum cut-off, A would be needed. Clearly 
2‘mT^ < vricr < A in order to allow fluctuations of the field 
on the scale of the pion compton wavelength. What is 
less obvious is that, with A of the order of ten or twenty, a 
cut-off much larger than 1 GeV leads to a negative effec¬ 
tive coupling, thus constraining 600MeV < A < IGeV. 
This is equivalent to a length- scale cut-off of the order of 
0.2fm and hence the model wouldn’t make much sense on 
time-scales much less than 1 fm / c. It is safe to assume 
that we should also not take the classical theory seriously 
on such length scales. 

At these sorts of energies, in the centre of mass frame, 
both the incoming nuclei appear highly Lorentz con¬ 
tracted into pancake shapes. In addition, experiments 
see uniform particle production as a function of rapidity, 
at least from the collision or central rapidity region. Both 
of these facts imply an approximate Lorentz boost invari¬ 
ance and consequently that physical quantities should 
depend only on the proper time r = y/t'^ — where x 


is the co-ordinate along the beam-pipe with zero at the 
collision point. The consequence of this is that initially 
at least the expansion of the plasma will be linear along 
the beam pipe. This should be true for times ( or dis¬ 
tances from the collision axis ) of the order of a nuclear 
radius, or t 1.2(Ai -|- « 7 fm / c for lead or 

uranium, where Ai and A2 are the mass numbers of the 
colliding nuclei. At later times, one would expect three 
dimensional rather than linear expansion but since this 
volume is at least as big as the region in which protons 
and neutrons are likely to be formed we ignore this. 

Since we are assuming, as is conventional, that the 
sigma model may be treated classically in the context 
of the chiral transition, it follows that during the initial 
linear expansion of the plasma, the temperature falls off 
like:- 


T = Ti 



a. 


where a depends on the speed of sound in the plasma and 
is 1/3 in the case of an ultra-relativistic plasma |^. The 
fact that the temperature can only depend on r implies 
that the plasma is hottest just behind the receeding pan¬ 
cakes and coolest in the central region, and cools from 
the inside outwards, somewhat like a baked-alaska. 


Important Time-Scales in the Plasma 



FIG. 2. A pictorial representation of the evolution of the 
plasma in space time. The remains of the incident heavy 
ions receed down the beam-pipe on light-like trajectories and 
all physical quantities and events are either specified or take 
place on space-like hypersurfaces of constant proper time. 

Glearly, the plasma and the sigma model can only be 
treated classically while there is more than one pion per 
pion compton wavelength:- 


kT > kTi 



rriT^c? 


or in other words up until the proper time th such that:- 
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Th < Tl 





For an initial temperature of about 200 MeV / fm^ this 
works out to be about 3Atj if one takes the pion Comp¬ 
ton wavelength to be one fermi and proportionally larger 
if one takes the value 1.5 fm. This is worth noting, since 
with the smaller value it is conceivable that the plasma 
would hadronize sufficiently early to be of dynamical in¬ 
terest. The relation between these time-scales is shown 
in figure ||. 


IV. ZUREK SCENARIO IN A QCD PLASMA 

To summarise, our picture of a heavy ion collision in 
the laboratory frame is of two highly lorentz contracted 
nuclei colliding to produce an approximately boost invari¬ 
ant plasma, initially in local thermal equilibrium. This 
then cools from the inside out, expanding initially lin¬ 
early, with correlation domains which grow as the phase 
transition is approached. This we will describe using a 
classical treatment of the linear sigma model, with criti¬ 
cal exponents calculated in the H = 0 limit. In order to 
predict the number of protons and neutrons produced in 
the central rapidity region, we need to know how many 
Skyrmion like configurations of the pion field will be pro¬ 
duced in traversing the phase transition. 

Although they will not actually be topologically stable 
in the sigma model, and will only become topologically 
stable objects when the plasma has cooled sufficiently be¬ 
low the transition that the field is almost always on the 
vacuum manifold and the Skyrme model is appropriate, 
the formation of these Skyrmion like field configurations 
will presumably be very similar to the formation of reg - 
ular topological defects by the Kibble mechanism [^ . 
In fact, counting the number of Skyrmions produced in 
this symmetry-breaking phase transition is equivalent to 
counting the number of topological defects produced in 
many other phase transitions, including those in super¬ 
fluid helium which have been used to experimentally test 
our ideas concerning defect formation. This is the cru¬ 
cial point which we exploit here in order to estimate the 
haryon density in the central rapidity region. 

In the present context, the Kibble mechanism would 
work as follows. During the transition, the 0(4) scalar 
pion field begins to fall from the false ground-state into 
the true ground-state, choosing a point on the ground- 
state manifold at each point in space. We will assume 
that the bias induced by the sigma term is small and 
that the point on the vacuum manifold is chosen approx¬ 
imately at random. Our sigma model has a second order 
phase transition so this collapse to the true ground-state 
will occur by phase separation, the resulting field config¬ 
uration being one of domains within each of which the 
scalar field has relaxed to a constant ground-state value. 


In the conventional Kibble mechanism, one then argues 
that continuity and single valuedness will sometimes force 
the field to remain in the false ground-state between some 
of the domains. This requires at least one zero of the field 
which would have topological stability and characterise 
a defect. The density of defects is then closely linked to 
the number of domains as shown in figure 



FIG. 3. In the conventional Kibble mechanism, topological 
defects, most often vortices of some kind, are formed between 
the domains of correlated field. 

In the case of the linear sigma model, however, it would 
be perfectly possible for the field to have Skyrmion like 
winding but still not be a defect since the field could 
unwind inside the surface over which the winding is com¬ 
puted. Counting protons and neutrons with the conven¬ 
tional topological current then might be hard. We can 
however still exploit heuristic evidence from simulations 
which indicates that there will be one real Skyrmion con¬ 
figuration formed for every twenty-five to a hundred do¬ 
main sizes [Q. Calculating the proton-neutron density 
then becomes a case of finding how many correlation do¬ 
mains are formed within the plasma. 

The question of how big the correlation domains are 
when defects are formed is the subject of the Zurek sce¬ 
nario for the formation of topological defects 0. This 
scenario provides an estimate of the initial defect den¬ 
sity immediately after the phase transition in a particu¬ 
lar case of the Kibble mechanism, namely a rapid quench 
through a second order phase transition. Although this 
scenario arose through considering the possible formation 
of cosmic strings in the early universe, it has since been 
tested using the superfluid transitions in both "‘He Q and 
^He [p|JlC|] and has so far been consistent with all obser¬ 
vations. In fact the systems in which the scenario has 
been tested are all strongly coupled in some sense and 
are hence far closer to the QCD transition in terms of 
their dynamics than they are to the sorts of cosmological 
phase transitions typically considered. 

In the Zurek scenario, the prediction of the defect den¬ 
sity formed during a symmetry breaking phase transi- 
tiondepends on the phenomenon of critical slowing down. 
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As the plasma cools and approaches the phase transi¬ 
tion, the correlation length grows and the relaxation rate 
of the pion field gradually decreases like r = At 

some stage during the cooling process therefore the re¬ 
laxation time-scale of the pion field will become longer 
than the time-scale on which the plasma is cooling. In 
other words, the pion field will no longer be in thermal 
equilibrium and will be unable to keep up with the cool¬ 
ing of the plasma]^ This is shown schematically in figure 

I- 


Correlation Length in the Zurek Scenario 



FIG. 4. Pictorial representation of the evolution of the cor¬ 
relation length according to the Zurek scenario. Initially the 
plasma is sufficiently close to equilibrium that the correla¬ 
tion length is able to track the quench quite well and the 
correlation length stays quite close to the equilibrium value. 
Eventually, as the transition is approached, however, it must 
go out of equilibrium at a particular proper time, tz- 

Hence, during the cooling of the QCD plasma, two 
regimes can be distinguished. Initially, sufficiently far 
away from the critical temperature the relaxation time- 
scale is much smaller than the time on which the cooling 
is proceeding and the pion field can maintain itself in 
local equilibrium with a correlation length ^(T). By con¬ 
trast, when the plasma has cooled to a relative tempera¬ 
ture in the vicinity of the phase transition, the pion field 
is effectively frozen compared to the time-scale on which 
the plasma is cooling. Thus, whatever the configuration 
and correlation length of the field at this relative temper¬ 
ature, it will be frozen in until after the phase transition 
when the field is again in thermal equilibrium. The time 
at which the plasma moves from one regime to another 
we will refer to as the Zurek time t. Quantities evaluated 
at this time will be denoted by a tilde. 


^Actually, the relaxation time of the field is mode-dependent 
so that, at least in a Gaussian approximation, different modes 
of the field would go out of equilibrium at different times. 
Remarkably, however, this turns out not to affect alter a naive 
estimation of the defect density 


Of course, by the time the field comes back into equilib¬ 
rium after the phase transition, the it may not be appro¬ 
priate to treat the pion field classically The important 
point, however, is that so long as the field behaves clas¬ 
sically until the Zurek time, when the field first goes out 
of equilibrium, the length scale which is imprinted on the 
plasma is the coherence length at the Zurek time. 

As an aside, we note that even if the chiral and de- 
confining transitions are actually separate and yet occur 
fairly close together, or even if there is not real transition 
but merely a slow change of phase, so long as there is a 
critical slowing down of the field relaxation time only one 
length will be frozen in. It may be possible to study the 
question of two transitions further in ^He or two dimen¬ 
sional networks of superconducting wire. 

Unlike cosmological phase transitions, however, which 
are always assumed to occur homogeneously, the quark 
gluon plasma cools from the inside out somewhat like a 
baked-alaska. In other words the second order transition 
into the symmetry broken phase occurs by the propaga¬ 
tion of a temperature front. The most striking example 
of this occurs in nematic liquid crystals where the phase 
front moves particularly slowly, but the transition pro¬ 
ceeds in the same way for many different systems. The 
speed of this front may be calculated as follows. 

Since the temperature depends on the proper time r = 
as T = at any given time, the 

temperature gradient will be:- 

^=TiaTf{cH^-x^)-^x 

ox 

Similarly, for any given position, the rate of change of 
temperature with respect to conventional time will be:- 

^ = —c^aT[Tf{c^t^ — x^)~^~ t 
dt ^ ^ ’ 

In other words, the length and time scales associated with 
temperature change will be:- 

1 1 dT ax 

At T dx 

1 1 dT c^at , , 

tq T dt 

Hence, in the centre of mass frame, this front will prop¬ 
agate with speed:- 


At c^t 



However, for the region inside the plasma, x < ct ex¬ 
cept very close to the receeding crepes where the exact 


Undeed, it will turn ont that the pion field goes ont of equi¬ 
librium only marginally before the end of the classical region. 
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boost invariance breaks down. Hence the acausal tem¬ 
perature front will always move at least as fast as the 
speed of light. In other words, the phase front propa¬ 
gates not only faster than the fluid which moves with a 
bulk velocity along the beam pipe oi xjt but also faster 
than the sound speed which is relevant for the equilibra¬ 
tion of the plasma. Effectively this means that there is 
no difference between this baked-alaska cooling and the 
homogeneous case as far as the formation of topological 
defects is concerned |^ . 

Let us therefore compute the relative temperature at 
which the textures will be frozen in for some point within 
the plasma with co-ordinates (x, t) in the centre of mass 
frame. Clearly, the field will go out of equilibrium when 
the time-scale associated with the cooling, tg, is equal to 
the relaxation rate of the plasma, Ir. As earlier, 

1 c^at 
tq 


and 


Suppose we consider a point in the plasma a fixed frac¬ 
tion \/l — of the distance from the point of impact to 
the position of the crepes at the outermost extremity 
of the plasma with x « ct. At any time t, this point 
(t, -^/l — f^ct) will have the proper time 


temperature everywhere and a single length-scale is im¬ 
printed on the plasma. 

The approximate Lorentz boost invariance of the 
plasma implied that physical quantities such as the initial 
conditions were specified on space-like Wpersurfaces of 
constant proper time, as shown in figure g. Although the 
plasma goes out of equilibrium at different times depend¬ 
ing on how close to the receeding crepes it is, the freezing 
in of topological defects actually occurs on a space-like 
hypersurface of constant proper time. 

Since for a; = 0 the proper time is equal to ct, us¬ 
ing the fact that T = T/(t//t)“ one finds that tc = 
Tj. Substituting our expression for temper¬ 
ature in terms of time at a; = 0 into our definition of 
relative temperature e = T/Tc — 1 one finds:- 


T 


Tl 


Tl 1 
Tc 1 + e 


l/a 


Rather than solving for e numerically, however, it is 
possible to obtain an approximate solution by specialising 
to the case of a; = 0 or equivalently / = 1 and using a 
slightly different criterion for when the pion field goes out 
of equilibrium, namely that the relaxation time is equal 
to the time remaining until the phase transition occurs:- 

tc - t = tR 


The criterion for going out of equilibrium then 
becomes:- 


T = x/c2t2 - a;2 ^ ^^2^2 _ ^2^2 ^ l ^ 

and consequently the quench time-scale at this point will 
be tg = ft/a. Since we know how to calculate the rel¬ 
ative temperature e as a function of proper time r we 
can relate ft to the relative temperature at the point 
{t, \/l — f^ct) in the plasma as:- 





-| l/c 


Equating the two time-scales then yields:- 


acto f Ti\ 


Ija 



This can be solved numerically for the relative temper¬ 
ature e which is frozen into the 0(4) field when it goes 
out of equilibrium and the length-scale imprinted on the 
plasma will be:- 



It will be seen that this relative temperature does not 
depend on /, the position in the plasma. In other words 
the plasma goes out of equilibrium at the same relative 



Exploiting the fact that e is likely to be small near the 
phase transition when the field actually goes out of equi¬ 
librium and in any case will certainly be smaller than 
the inital value of roughly 0.25 we can series expand 
(1 -I- to give:- 


~l 

h 


1 

1 + vz 

6|x=0 — 

a — 
tl 

J 



Since we know that the plasma freezes out at the same 
relative temperature everywhere, and in particular it 
freezes out on a specific space-like hypersurface of con¬ 
stant proper time, we may as well use this approximate 
solution rather than solving the previous equation nu¬ 
merically. 

Using this approximate solution for the relative tem¬ 
perature e, we can calculate the length scale which is 
frozen into the plasma during the phase transition, the 
Zurek length. Since we know that volume of the QCD 
plasma at any particular time we can then calculate how 
many correlation volumes are frozen into the plasma ac¬ 
cording to the Zurek scenario and hence how many pro¬ 
tons and neutrons we would expect to see in the central 
rapidity region. 



















Since the initial conditions are not well known the re¬ 
sults for a variety of plausible initial temperatures and 
proper times are shown:- 


Ti/fm 

Ti / MeV/fm^ 
etc / fm 
e/ 

1 

200 

2.0 

0.25 

1 

300 

6.6 

0.88 

2 

250 

7.6 

0.56 

4 

200 

7.8 

0.25 

4 

300 

26.4 

0.88 

e 

1 / fm 

Plasma Volume / fm^ 
Number of Domains 

0.36 

1.5 

246 

18 

0.20 

2.3 

1170 

24 

0.19 

2.4 

1386 

25 

0.19 

2.4 

1417 

26 

0.10 

3.8 

6098 

28 


where we have taken the QCD plasma to be a cyliner 
of width 14fm and length equal to 2ci. Although both 
the Zurek proper time, or equivalently the frozen in do¬ 
main size, and the volume of the plasma at this time 
are quite sensitive to the initial conditions, which are 
not well known, the total number of correlation volumes 
frozen into the plasma is much less sensitive since it is 
the ratio V /tt^^ . 

In other words there will be somewhere between 15 and 
30 coherence volumes frozen into the plasma. From sim¬ 
ulations, however, we expect to get roughly one skyrmion 
configuration per every 25 - 100 domains. In order to 
create a proton-antiproton pair then ^ we would need to 
have between 50 and 200 coherence volumes and also be 
lucky enough that the pair didn’t immediately annihilate. 
Avoiding annihilation is not entirely implausible since the 
distance between the nucleon and anti-nucleon must be 
at least as great as the Zurek length which in this case 
is slightly greater than the proton Compton wavelength. 
Forming a sufficiently large number of domains that the 
resulting proton and neutron density would be observ¬ 
able on the other hand is far more difficult. Roughly 
speaking, however, this result indicates that the central 
rapidity region of the plasma is likely to remain free of 
protons and neutrons as is usually thought. 

With the smallest possible length-scale imprinted, 
namely the pion compton wavelength, the situation is 
somewhat better and one could conceivably see a rea¬ 
sonable number of protons and neutrons. However, this 
corresponds to almost no domain growth and in this con¬ 
text DCCs would also be ruled out. 

It is possible, however, that we might be extraordi¬ 
narily lucky since the predicted values are necessarily 
rather approximate and do not categorically exclude the 
possibility of a detectable baryon number in the central 
rapidity region. Also, the freeze out occurs at a suffi¬ 
ciently large proper time that it is debatable whether we 
are still justified in treating the model classically. The 
model is already known to be flawed since it does not 
conserve baryon number and the results are somewhat 


^Creation of a single proton would violate baryon number 
conservation. 


suspect because of this. It is certainly possible to do a 
slightly better calculation to include the quantum me¬ 
chanical aspects of the theory, but this would not cure 
the problem of baryon number conservation. It seems 
likely that the only way to make a significant improve¬ 
ment would be to improve the model somehow, and this 
would probably make further analytic work intractable 
and necessitate a simulation. 


V. CONCLUSIONS 

The conventional wisdom is that the central rapidity 
region of a heavy-ion collision in which a quark-gluon 
plasma is produced should be almost baryon free. A di¬ 
rect estimate of this initial baryon density in the context 
of the 0(4) linear sigma model shows this belief to be 
well founded and entirely consistent with the production 
of dis-oriented chiral condensates. 

However, a non-zero baryon density is not ruled out 
by many orders of magnitude and the result is some¬ 
what sensitive to the initial conditions. It is therefore 
conceivable that with a slightly larger volume of QCD 
plasma one might occasionally produce a small but non¬ 
zero baryon density. To determine whether it is worth¬ 
while looking for this experimentally further theoretical 
work is necessary. 
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